On the Equivalence between Mode-Coupling Theory and Replica Approach for 

Super-Cooled Liquids 
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We show that the replica approach to glassy dynamics provides, in spite of its static nature, a 
characterization of critical dynamics in the /3-regime of super-cooled liquids that is equivalent to the 
one of Mode-Coupling-Theory, both qualitatively and quantitatively. The nature and extent of this 
equivalence is discussed in connection to the main open problems of the current theory. 



It is nowadays well established that the Modc- 
Coupling-Theory (MCT) yields and accurate description 
of the early stages of the dynamical slowing down in 
super-cooled glass- forming liquids One of its main 
successes is the prediction that the slow relaxation has a 
two-step nature with time correlators developing plateaus 
that are expected to diverge at some critical temperature 
characterized by a transition to a glass phase. 

The main drawback of the theory is that the crit- 
ical temperature is actually significantly larger than 
the one of the glass-transition in numerical simulations 
and it seems rather to mark a dynamical cross-over 
Nevertheless various quantities, like the so-called 
non-ergodicity parameter (essentially the height of the 
plateau of the correlators), compare very well with the 
numerical data. Furthermore the theory provides a de- 
tailed set of predictions for the behavior of the critical 
correlators and the dynamical exponents Q that repro- 
duce quite well the numerical data (j-|6(. For these rea- 
sons many believe that the inclusion of some sort of 
corrections to the theory, possibly taking into account 
finite-dimensional effects leading to some some sort of 
hopping processes, could extend MCT down to the lab- 
oratory glass transition. Unfortunately due to the un- 
controlled nature of Sjogren's approximation Q which is 
at the heart of standard MCT [8| it is not clear how to 
compute systematically corrections to it. 

On the other hand it has been discovered by Kirk- 
patrick, Thirumalai and Wolynes [9( that a certain class 
of mean-fied spin-glass models exhibits the same dynam- 
ical features of MCT. This result is important for various 
reasons. It shows that MCT works for systems that are 
completely different at the microscopic level from a liq- 
uid hinting that it might have some universal features. 
Furthermore it allows to include, at least phenomeno- 
logically, finite-size effects in the theory. Indeed crit- 
ical dynamics in these models can be tracked back to 
the splitting of the equilibrium state into an exponential 
number of metastable states and the use of nucleation 
arguments on this mean-field picture has led to the for- 
mulation of the so-called Random-First-Order-Theory of 
the glass transition [H, [ljj. Another important conse- 
quence is that some features of critical dynamics of these 
models, namely the critical temperature and the non- 
ergodicity parameter, can be computed without solving 
explicitly the dynamics by means of an essentially a static 



approach: the replica method. This observation has led 
to the application of the replica method directly in the 
context of liquid theory, leading to a novel set of quanti- 
tative predictions for the non-ergodicity parameter and 
for the critical temperature [11} . 

In the replica approach one is typically concerned with 
a replicated order parameter q a b that describes the ap- 
propriate microscopic two-point correlation {i.e. density- 
density in liquids or spin-spin in spin-glasses) between 
two identical replicas of the system. Using standard 
statistical mechanics procedures one derives a (usually 
approximate) equation for q a ^ and the dynamical glass 
transition is characterized by the abrupt appearance of 
a novel solution of the equation characterized by a value 
of q a b different from the value corresponding to the high 
temperature phase (say the liquid or paramagnetic phase 
depending on the system considered). The solution de- 
pends on the number m of replicas considered and for 
technical reason one should take the limit m — > 1 in 
order to obtain the correct result. In the structural 
glass context it is customary to use methods based essen- 
tially on the Hyper-netted-chain (HNC) approximation 
[ill fl2| that give predictions for the non ergodicity pa- 
rameter f(q) quantitatively different from those of MCT. 
Because of this discrepancy it has been assumed in the 
past that the two methods are intrinsically different. In- 
stead Szamel [l3| has shown that starting from replicated 
Ornstein-Zernicke equation and closing them with a suit- 
able approximation scheme one can recover exactly the 
results of MCT. Therefore it can be argued that the two 
approaches predict the same physics and that discrepan- 
cies arise only when inequivalent approximation schemes 
are used. 

Obviously the two methods cannot be completely 
equivalent. In principle MCT can describes dynamics on 
both small and large times scale while the replica method 
is intrinsically a static approach. Therefore it is reason- 
able that the latter can describe a static (or rather quasi- 
static) quantity (like the ergodicity-breaking parameter) 
while it cannot obviously tell anything on short-time dy- 
namics. Previously it was assumed that it could not tell 
anythin g on dynamics at all. Instead it has been recently 
realized [14j that it can be used to characterize the criti- 
cal behavior on the large time scale of the (3 regime, i.e. 
the plateau in the time correlators. In order to do so one 
has to perform an expansion of the equation for the order 
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parameter q ab around the solution q ab — q c corresponding 
to the glassy phase. Independently of the approximate 
method used in order to obtain the equation it follows the 
we should get an equation in powers of q a b = <Zc + 5q ab of 
this form: 

= T + m 2 Sq bc + 5q ac J + m 3 ^ &<lcd + 
V c / cd 



+ w 1 (5q 2 ) ab + w 2 Sq 2 ab + 0(dq 3 ) 



(1) 



The various coefficients depend on the given model and 
on the external parameters, in particular t vanishes lin- 
early with the external parameters at the critical point 
and is positive(negative) in the glassy(liquid) phase. In 
[lij it has been shown that the structure of the above 
equation determines the behavior of the time correlators 
near the critical point, in particular the critical behavior 
of the two-point time correlation C (t) around the plateau 
value C p = q c is given by: 



5C{t) EE C(t) 



C P = M 1/2 /±(*/**) t*oc 
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where the function /_ has to be chosen in the liquid 
phase (r < 0) while the function / + has to be choosen in 
the glassy phase (r > 0) and t* is the time scale of the 
/3-regime. The function f±{x) obeys the scale-invariant 
equation: 

±1 = fl{x) (l - — ) + l\f±(x-y)-f±(x))f ± {y)dy 

(3) 

For small values of x both the functions f± (x) diverge as 
l/x a , while for large values of x f+(x) goes to a constant 
while f-(x) diverges as — x b where the exponents a and 
b are determined by the so-called parameter exponent A 
which turns out to be equal to the ratio wijw\\ 



Aee^ = 

Wi 



r 2 (i-a) _ r 2 (i + &) 

T(l - 2a) ~ r(l + 26) 



(4) 



Note the the above results are quite general and the mi- 
croscopic details of the given model enters only in the 
actual value of the parameter exponent A. The above 
equations are equivalent to those obtained originally in 
the context of MCT by Gotze (see eqs. 22a and 24 Q) 
who also recognized their exact nature, i.e. the fact that 
the (model-dependent) value of A is the only quantity 
sensitive to the approximation scheme used. Note the 
this result is another hint at the equivalence between the 
replica approach and MCT because it shows that their 
predictions do coincide when they are exact. 

The result derived in |14j provides a simple recipe to 
determine A within the replica approach: one has to ob- 
tain a closed equation for the non-ergodicity parameter 
using his preferred approximation scheme, expand the 
equation around the critical value of q a b up to second or- 
der as in eq. ([!]), read the coefficients w 2 , w\ and finally 



identify A with their ratio. This recipe has been applied 
to a number of mean-field spin-glass models yielding new 
analytical predictions concerning critical dynamics [TEl - 
[Tsl ] . When there was already an explicit solution of the 
dynamics the novel computations offered an a posteri- 
ori validation of the recipe, while in other cases the new 
predictions were compared with existing numerical sim- 
ulations. More recently the recipe has been used in the 
context of the replicated theory of super-cooled liquids 
based on the HNC approximation [12| and for a wide 
range of models a (quite reasonable) value of A « .7 was 
obtained. In order to show that the replica approach 
and MCT provide equivalent descriptions of the long- 
time critical dynamics in the following the recipe will be 
instead applied to Szamel's approximation scheme. 

The issue of the equivalence between the two ap- 
proaches is not an academic one. On one hand MCT 
predictions are more accurate at present than those ob- 
tained with the HNC approximation, at least in three di- 
mension. On the other hand there are situations in which 
the replica method may have considerable advantages. 
Indeed in order to derive many quantities within MCT 
one takes the large-time limit in the dynamical equations. 
The replica approach can be viewed as a compact way to 
describe the algebra of this limit where different replicas 
correspond to configurations visited dynamically at dis- 
tant times. Certainly the main open problem concerning 
MCT is the systematic inclusion of corrections and the 
role finite-dimensional effects/finite-size effects 
The replica approach offers the possibility to study this 
phenomena in a field-theoretic framework that is defini- 
tively simpler than the dynamical one. Surprisingly it 
has been recently discovered that the loop corrections to 
this theory can be controlled at all orders and a series 
of novel unexpected predictions have been obtained [l9| . 
However this study confirmed that the phase transition 
predicted by MCT is only a dynamical cross-over and 
therefore the purely static replica approach, as expected, 
is inconsistent beyond perturbation theory. The theory 
needs therefore to be supplemented with some form of 
critical dynamics but unfortunately in spite of impor- 
tant progresses [2lJ we are still lacking full control of the 
critical behavior of the dynamical propagators, like the 
well-known \<t susceptibility. 

Computation of the parameter exponent Our start- 
ing point is the replicated Ornstcin-Zernicke equations 
introduced by Szamel in Ref. [13| . 

Sn ab (k) = n 2 c ab (k) + n^2c ac (k)Sn cb (k) (5) 



where 5n ab (\i) is the Fourier transform of Sn ab (ri, = 
"ab(ri,r2) — n 2 , i.e. the non trivial part of the two- 
particle density, n being the total particle density. The 
above equations can be closed provided an expression 
for the direct correlations c ab (k) in terms of the 6n ab (k) 
itself is given. The single replica correlation functions are 
associated by the definition to the static structure factor 



3 



according to: 



Sn aa (k) = n(S(k) - 1) 



(6) 



Under the assumption of Replica-Symmetry (RS) the 
equation for the diagonal components factorizes from 
that of the off-diagonal components in the limit m — > 1 
and equation ([5]) reduces to 



Sn aa (k) = n 2 c aa {k)S(k) 



(7) 



the above equation combined with the definition ([6]) can 
be seen as a definition of the direct correlation function: 



,(k) = c(k) 



S(k) - 1 
nS(k) 



(8) 



Indeed the static structure factor is usually considered 
an input in MCT computations. Szamel has shown that 
that making an appropriate set of approximations in the 
context of the Yvon-Born-Green hierarchy one obtains an 
expression for the off-diagonal direct correlation identical 
to Sjogren's vertex Q. More precisely if we define the 
replicated non-ergodicity parameter / a b(k) according to: 



5n ab (k) =nf ab (k)S(k) 



(9) 



we obtain the expression: 



Cab(k) = J dqdpG(k,q,p)f ab (q)f ab (p) (10) 

where the function G(k, q, p) has precisely Sjogren's 
form: 



G(k,q, p) 



2k 2 (2tt 



-S[k-p- q }S(q)S(p)(k-[qc(q)+pc(py_ 

(11) 

Now we multiply cq. ([5]) by the inverse of the matrix 
nl + Sn ac (k), (I is the identity in replica and momentum 
space), and we use the above equation to obtain: 

r 5n , ) (k)=n [ dqdpG(k,q,p)f ab (q)f ab (p). 

In order to perform an expansion around the RS solution 
we rewrite: 

Sn ab (k) = 5n* b s (k) + nS(k)(l - f(k)) 2 Sq ab (k) (13) 
where 5n^ b s (k) has a RS structure 



5<f(k) = n(S(k)-l) 
6n% b S (k) = nf(q)S(q) a^b 



(14) 
(15) 



The factor nS(k)(l — f(k)) 2 has been used in the defini- 
tion in order to make contact with ref . Q . According to 
[l4j we have to expand equation (TP2)) in powers of Sq ab (k) 
assuming Sq aa (k) = 0. 

In order to perform the expansion we first note that 
the r.h.s. (where the crucial approximations are made) 



has a simple diagonal structure in replica space and only 
the l.h.s. has a non-trivial structure. The zero-th or- 
der term can be computed using standard properties of 
replica symmetric matrices and leads after some algebra 
to: 



5 n 



nl + 5n 



(k) 



/(*) 



S(k)(l - f(k)) 



Q{5q), a^b 



(16) 

In order to compute the next terms we note that in 
a power series expansion they reduce to the computa- 
tion of objects of three possible forms: 1) (SnRs)^ b , 
2 ) ( Sn RS S 1 Sn Rs)ab and 3) (Sn r RS Sq Sn s RS Sq Sn RS ) ab for 
general integer r, s, t (all products and powers are inten- 
dend in matrix sense and we have not written explic- 
itly the replica and Fourier indexes). At linear order 
the expansion generates terms of three types: miSq ab , 
TO 2Sc(^9ac + Upland m^^cd^lcd- However, as ex- 
plained in Ref. [1J|, only the first one, that depends 
simultaneously on both indexes a and 6, is relevant 
for the dynamics [23| . The computation of this term 
is simpler than the others because since both indexes 
must be present every time we have a multiplication by 
Sn RS (k) = n(S(k)(l - f(k)) - l)5 ab + f(k)S(k) we are 
forced to take the term with the Kronecker delta. The 
same observation is correct also at quadratic order in 
5q(k), i.e. for the term ^ c 5q ac (k)5q cb (k). Therefore 
in order to compute the terms relevant for the dynam- 
ics we can replace the much simpler form Sn ab (k) = 
n(S(k)(l - f(k)) - l)S ab + nS(q)(l - f(q)) 2 8q ab (k) into 
the l.h.s. ofeq. (fT^j) . Considering only the terms relevant 
for the dynamics, the final expression is: 



^,2 



/(*) 



Sq ab (k) + (1 - /(fc))(^ 2 ) ah (k) + 



(WW) 
+ nS(k) J dqdpG(k,q,p)/(g)/(p) 



2nS(k) / dqd P G(k,q,p)(l-/(< Z )) 2 ( 5 (?ab (q)/(p) 



+ nS(k) J dqdp G(k, q, p) x 

X (1 - f(q)) 2 {\ - f(p)) 2 6q ab (<l)5q ab (p) 



(17) 



The equation at zero-th order determines the non- 
ergodicity parameter while the next orders determine the 
critical temperature and the parameter exponent. The 
critical temperature is fixed by the condition that there 
is a direction for 5q ab {k) such that the coefficient of the 
linear term vanishes. This corresponds to the condition 
that the following linear operator has a zero eigenvalue 
H = 0: 

= -tf(k) + 2nS(k) J dqd P G(k,q,p) X 
x (l-f(q)) 2 tf(q)f(p) (18) 

where £/?(^) is the right eigenvector corresponding to 
eigenvalue /i. Note that this condition is equivalent to 
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the condition that the derivative with respect to f(q) 
of the zero-th order term vanishes and therefore marks 
the point in temperature where a non vanishing f(q) ap- 
pears. This equivalence between two different conditions 
is a consequence of the m — > 1 limit and of the fact that 
in this limit the so-called replicon and longitudinal eigen- 
value are degenerate leading to many non trivial features 
concerning fluctuations [I'.i . The behavior at the critical 
point is controlled by the critical eigenvector £,J}(k) and 
we may write at leading order: 



Note that the above definition must be supplemented 
with a normalization condition on the eigenvector, in or- 
der to make contact with ref. Q we adopt the following 
normalization for the right and left eigenvectors: 



dk$(k)tf(k) = i, dktf(k)(i-f(km*(k)$(k) = i 



Sq a b(k) = 8q ab £,o(k) 



(19) 



Multiplying the quadratic part of equation (|17|) by the 
left eigenvector £q (k) and integrating we obtain the fol- 
lowing expression: 



J dk^(k)(l 



(Sq ab ) 2 



I 



dkdqdp#(*)S(fc)G(k J q,p)(l - f(q)) 2 (l ~ /(p)) 2 ^(<z)C(p) 



m 



Note that the first coefficient is fixed to one by the nor- 
malization. According to [Tij the parameter exponent is 
given by the ratio between the two quadratic coefficients 
leading to: 



A = 



dkdqdp Co (k)S(k)G(k, q, p) 



X (l-f( q )) 2 (l-f(p)) 2 tf(q)tf(p) 



(21) 



Identifying the two-point vertex in [3] as 

2nS(k) G(k, q, p) we see that that eqs. (fT8|) and (|21|) arc 

respectively equivalent to eq. (13a) and (15a) in [3j. 

Conclusions The recipe of |14j has been applied to 
the approximation scheme proposed by Szamel in order 
to show that the replica approach gives a characteriza- 
tion of the /3-regime of super-cooled liquids that is equiv- 
alent to the one of MCT both qualitatively and quanti- 



tatively provided the same approximations are used. A 
novel result obtained within the replica approach is that 
an equation like ([1} follows from a replicated free energy 
14]. This allows to connect the coefficients wi and W2 
with two six-point susceptibilities lj\ and U2 which are 
in principle measurable in a quasi-static framework. Ac- 
cording to [lH their ratio is equal to A and therefore eq. 
(|2ip fixes their value up to an unknown constant. The 
computation of this constant is left for future work. In 
the present framework it could be determined consider- 
ing the response of the non-ergodicity parameter to a 
random pinning field along the lines of [22j. 
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